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Abstract 

We use the r-matrix formulation to show the integrability of geodesic 
flow on an iV-dimensional space with coordinates qt, with k = 1, N, 
equipped with the co-metric <? lJ = e~' Qi ~ *l(2 — e~' 9i ~' 2j '). This flow is 
generated by a symmetry of the integrable partial differential equation 
(pde) mt + um x + 3mu x = 0, m = u — a 2 u xx (a is a constant). This 
equation - called the Degasperis-Procesi (DP) equation - was recently 
proven to be completely integrable and possess peakon solutions by 
Degasperis, Holm and Hone 0. The isospectral eigenvalue problem 
associated with the integrable DP equation is used to find a new L- 
matrix, called the Lax matrix, for the geodesic dynamical flow. By 
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employing this Lax matrix we obtain the r-matrix for the integrablc 
geodesic flow. 
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1 Introduction 

The 6 = 3 peakon equation and its isospectral problem 

We begin with the case b = 3 of the b— weighted peakon equation. This is the 
evolutionary equation defined on the real line as, 

m t + um x + b mu x = , m = u — a 2 u xx , lim m = 0, (1-1) 

in which the subscripts denote partial derivatives with respect to the indepen- 
dent variables x and t. For any values of the dimensionless constant b > and 
constant lengthscale a, this equation admits exact iV-peakon solutions 

N 
i=l 

in which the 2N time-dependent functions Pi(t) and i = 1,2, . . . , N, 

satisfy a system of ordinary differential equations whose character depends on 
the value of the bifurcation parameter b. The case b = 2 is the dispersionless 
limit of the integrable Camassa-Holm (CH) equation that was discovered for 
shallow water waves in H . As shown in Q , the CH equation with dispersion is 
one full order more accurate in asymptotic approximation beyond Korteweg- 
de Vries (KdV) for shallow water waves, yet it still preserves KdV's soliton 
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properties such as complete integrability via the inverse scattering transform 
(1ST) method. 

An equation equivalent to the case 6 = 3 of the peakon equation 
was singled out for special attention among a family of related equations by 
Degasperis and Procesi in 0. The peakon equation ( [Lip was shown to be 
completely integrable for the case 6 = 3 by Degasperis, Holm and Hone M, 
who found the following Lax pair consisting of a third order eigenvalue problem 
and a second-order evolutionary equation for the eigenf unction, 

ipxxx = -^ip x - Xrnip, (1.3) 
or 

1 2 

i[>t = ^ripxx-uipx + {u x + (1.4) 

crA 3A 

Compatibility ip xxx t = iptxxx implies Eq. (|1 . 1| ) with 6 = 3 provided dX/dt = 0. 

Thus, equation (|1.1| ) with 6 = 3, namely, 

m t + um x + 3mu x = , m = u — a 2 u xx , (1-5) 

was shown to be integrable by the inverse spectral transform for the isospectral 
eigenvalue problem ( |1.3|) and to possess an infinite sequence of conservation 
laws in Degasperis, Holm and Hone . Henceforth, we shall consider the case 
for a = 1 for the sake of simplicity in writing. (The lengthscale a is easily 
restored from dimensional considerations.) 

The first few conservation laws are given in the notation of M with a = 1 

by 

H-i = \ f u 3 dx , Hn = fmdx, . 

6 J ' (1.6) 

H x = y(v 2 xx + 5v 2 x + Av 2 )dx, H 5 = Jm^dx. 

We shall pay special attention to the quadratic conservation law Hi, in which 
the quantity v is defined as 

v := (4 - dlY l u = (4 - " dl)~ l m . (1.7) 
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Lax matrix for iV-peakon dynamics 

Substituting the iV-peakon solution, 

N N 

u(x,t)=J2P^ e ~ lX ~ qm > rn{x,t) = 2Y,Pj{t)^-^)\ (1-8) 

into the isospectral eigenvalue problem ( |1.3[ ) yields P 
1 1 " 

-^(x, t) = -^[l + ^(z - q 3 {t)) (1 - e -l*-«(*)l)] p^( qj (t)) . (1.9) 
i=i 

Setting ip(qi(t),t) = ipi{t) then gives the following matrix eigenvalue problem, 

2 N 

—*Pi = ^2 L ij^j, (1-10) 

where 

4- = [1 + s^nfe - <&■) (1 - e-l*-*')] Pj . (1.11) 

Let L denote the N x N matrix Ly. In Ref. 0, the authors used the two 
conserved quantities trL and trL 2 to solve the 2-peakon subdynamics of the 
the iV-peakon dynamics qk,Pk, with k = 1, ...,N, a = 1, satisfying 

Pj = ^T,k=iP3PkSgn{q j -q k )e'^-^\, 
Qi = Eti^e- 1 *-* 1 - 

Amongst other results, the authors in || discovered the two-peakon collision 
rules for N = 2 and gave explicit formulas for its phase shifts as functions of 
the asymptotic speeds of the two peakons. 
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2 A geodesic flow for 6 = 3 peakons 



The quantity used for determining the two-peakon N = 2 collision laws in ||, 
1 1 N 

H x = -trL 2 = -£ /'</',' " ''' ( 2 - e-1*-«1) , (2.1) 
is also the quadratic conservation law Hi in ( |1.6|) for the 6 = 3 peakon equation 



( |1.5|) , when Hi is evaluated on the N— peakon solution ( |1.8| ). 

The canonical Hamiltonian dynamics generated by Hi is geodesic motion 
on an iV-dimensional space with co-metric = e _ ' 9< ~ 9j '(2 — e"' 9 * - *'). As we 
shall show by finding its r— matrix structure in the remainder of the present 
paper, this geodesic motion canonically generated by the conservation law 
Hi = TrL 2 in Q2.1|) provides a new 2iV-dimensional integrable system, 

BH N 

» = If 1 = HPi e ~ Wk ~ qA ( 2 - c-'*-* 1 ) , (2.2) 
BH N 

Pk = -^ = -2p k ^p ] sgn{q ] -q k )e~\^\{l-e-^-^ . (2.3) 

We emphasize that these geodesic ifi-dynamics for Pk,Qk, are not the same 
as the iV— peakon dynamics in ( |1.12| ). Rather, we are studying the canonical 



Hamiltonian dynamics on the restriction to the peakon sector of the Hi -flow in 
the hierarchy of integrable equations associated with the isospectral problem 
for equation (|1.5|). 



i?-matrix results for the geodesic i^i-dynamics 

To show integrability of the canonical dynamics in (|2.2|) and (|2.3|) produced by 



Hamiltonian fl2.1| ), we shall use the r-matrix method, which reveals that trU 
for a certain Lax matrix L are constants in involution IHl . 
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To find the r-matrix structure for these ifx-dynamics for Pk,Qk, we shall 
introduce an alternative Lax matrix for the peakon dynamics of Eqn. ( |1.5|) , 



N 



I J ■ 



(2.4) 



where 



VPiPj A 



A(qi - qj) = yj(2 - e-\^-^\)e-\^-^\. 



(2.5) 
(2.6) 



The Lax matrix (12. 41) also satisfies 



1 1 N 

H t = -trL 2 = ^J2 PiPj^ qi ~ qA ( 2 - e-'*-* 1 ) 



(2.7) 



which is the Hamiltonian for the canonical dynamics in Eqs. ( |2.2|) and (|2.3|) . 
In Eq. (p.6|), we have 



and the function A(x) satisfies the following relations, 

1 - 

A(x) = -sgn(x) 2 _ |x| A(g), 
I - 4 . . 4 . . _ i 

4,- = - fc) = -A'fo - ft ) = -A' 5i , A' u = 

(-^ + ^)A(x)A(y) = A'(x)A(y)+A(x)A'(y) 



= -A{x)A{y) 

d d 
ox ay 



1 - e~ |x '! 1 - e~\ y \ 

sgn(x) 2 _ e _ ]x] + sgn(y) 2 _ e _ ]yl 



0. 



(2.1 



(2.9) 
(2.10) 
(2.11) 



(2.12) 
(2.13) 



We shall work in the canonical matrix basis E+j, 



(Eij)ki = SikSji, i,j,k,l = 1,...,N. 

To find the r-matrix structure for the ifi-dynamics in Eqs. ( |2.2| ) and ( |2.3| ), we 
consider the so-called fundamental Poisson bracket [[?]: 



{L U L 2 } = [r 12 ,L x ]-[r 2X ,L 2 ], (2.14) 



where 



In 



N 

E 

8, t 

N 



N 




1 = LijEij 


(8) 1, 






N 








k,l=l 




1"ij;stEij ® E st , 





JV 

{-^1,-^2} = {Lij, L k i}Eij ® E k [. 

i,j,k,l=l 

Here {Ly, L/^} is the standard Poisson bracket of two functions, 1 is the N xN 
unit matrix, and the quantities rij- st are to be determined. In Eq. ( |2.14j) , [•, •] 
denotes the usual commutator of matrices. 



After a lengthy calculation for both sides of Eq. ( 2.14 ), we obtain the 



following key equalities (whose detailed verification is given in the Appendix): 
[ r ih L \ii = 0, 
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[ r 3h L \lj 

[ r ihL\ij 
VihL} kj 

[ r jh L \lk 
[ r 3h L \kl 



[ r l3i L \lj = 0. 

[ru,L]ji = -y/PW A 'ji> 

Vih L\u = -y/PjPiA'ji, 

1 



0, for different s,t,j,l. 



^VPkPjiAjiAik)' , j^l,k;k^ I, 



where r j7 = Y l k, m r km,ji E km, r u = Y.k, m r km,uE km , are two N x N matrices 
whose entries are to be determined, L is the Lax matrix, and [•, L] fci stands for 
the fc-th row and the l-th colum element of [-,!/]. 

In matrix notation, all the above equalities can be rewritten as 



[r jh L] = &\ j^l, 
[ru,L] = B u , 

where B^ 1 , B a are the following two N x N matrices: 

/ ... l^/pmiAijA^) 1 ... l^pTpJiA.iAnY 



(2.15) 
(2.16) 



B jl 



\^pm{Ai A 0l )' 



\^/mrMjiA n y 



-VPJPiA'ij 



-VpJpi a '] 



.a 



\^pWp]( A 31 A in)' 



\y/PNPl{AijA jN )' ... \^/PnT]{ A 31 A In)' 
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and 



/ 



N! 



V 



-y/PNPlA'j 



N! 



J 



By solving Eqs. ( p,15| ) and ( |2.16| ), we have the following r-matrix structure: 



A' 



3,1=1 V « 



1 N 




-Ay 

+ Eij) + -r-Eu 
A lj 



A' k3 A kl (A^y 



E 



j.i 



AkjAij 



A 



ij 



Eu ® E 



;2.i7) 



Perhaps not unexpectedly, this non-contant r-matrix for the geodesic Hi- 
dynamics differs from the constant r-matrix associated with the CH equa- 
tion (b = 2) discovered by Ragnisco and Bruschi in M. According to the 
r-matrix method [Q], the auxiliary matrix M corresponding to L may be also 
constructed. But, in our case the r-matrix is dynamical, and therefore the 
construction of M is a rather complicated process, which we defer to another 
time. 



Concluding remarks 

In this paper, we found the r-matrix formulation for the integrable geodesic 
motion generated canonically by the quadratic quantity Hi in ( |2.1| ). This 
quantity arises by restriction to the peakon sector of a quadratic conservation 
law in the hierarchy of integrable equations associated with the isospectral 
problem for equation ( |1.5|) . The quadratic quantity Hi is also conserved for 
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the 2-peakon dynamics of the 1+1 integrable partial differential equation (|1.5| ) 
that was singled out in || and was proven to be completely integrable by 
the isospectral transform in |J. We also introduced a new Lax matrix L for 
the iV-peakon flows of the integrable equation ( |1.5| ) that facilitated the r- 
matrix calculations and for which Hi = |trL 2 . In later work, we shall discuss 
additional flows in the hierarchy of integrable equations associated with the 
isospectral problem for equation ( |1.5| ) and study their relationships to classical 
finite-dimensional integrable systems H. 
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Appendix 



The following computations are needed in verifying Eq. ( |2.14| ). 
First, we calculate the left hand side of Eq. ( |2.14| ). 



dq m 
dL kl 

dp m 



^/PiPjA'ijiSim - 5 jr 



A 



ki 



N , 

{Lij,L M } = / A 

m — -\ \ 



dLij dLki 



dLki dLij 



J V dq m dp m dq m dp n 



1 N 



m=l 



A, 



ki 



(Sim - $jm)(Pl$km + PAm) 



~ VPk~PlA' h 



A, 



ki 



(hm - 8i m )(pj5 im +Pi6j 



■jm) 



- Vpipi{AijAji)'S k j + y/PkPj{AijA ki )'ki\ 
\ WwPi{AkjAji)'8ji + ^/ptp~j(AiiA i: j)'5 ik 
+ y/PiPiiAijAjiYSkj + y/pkPj(A k iAij)'8ii\ , 



where the superscript / means Eq. ( 2.12j ) with the argument. 
Thus, we obtain the following formula, 



N 



{L±, L 2 } = 2_j {Lij,Lki}Eij (g) E k i 

i,j,k,l=l 
1 - 

= 9 $Z [y/PkPjiAkiAijYEjt ® E H + y/pkp^AkiAijYEij ® E, 



j,k,i=i 



12 



+ y/PkPj(A k iAij)' Ejt <g> E lk + y/PkPj{A k iAij)'Eij ® E kl ] 
1 - 

n ^ v/PfcWO 4 *^;)' ® £*j + #y ® Si* + E iX ® E lk + Eij ® £, 



iV 



j,k,l=l,j^k,l;k=£l 
N 

fc,Z=l 



(2.18) 



Next, we compute the right hand side of Eq. (p. 14 ) 
[ri 2 ,Li] - [r 21 ,L 2 ] 



A? 



X r^L*, [(Eij ® £ st ) (£ w <g> 1) - (£ fci ® 1) (E y ® £ sty 

i ,j,s,t,k,l=l 

-(E st ® £y)(l ® + (1 ® ® 
^ [r i:) - ;sfc (Lji(Eu <g> i£ s &) — Lu(Eij <8> £? s fc)) — ry ;ss L fc2 <S> EijEki) — (E ss (g> E k iEij))] 

i,j,s,k,l=l 
N 

^ [ry; S fc (Lji(Eu <g> £ sfc ) - L M (E y <g> £ sfc ))] 

N 

+ X h r j;;ss-M-£ss <8> + rij. ss L u (E ss <g> £y)] 
^ [rji-stLi^Eji <g> E st ) - Tij-gtU^Eij <g> £ st )] 
+ X [rji. ss L u (Eji <g> £7 SS - £ ss <g> - r ij]SS Lii(Eij <g> £ ss - £ ss ® £? y ) 

i,j,s,l,Sitj,l,j^l 

;ljLi s Ti s -ijL s i )E SS ® Ey)] 

+ X^ [( r u-ijLii — ru-ijLn + Tij-uLu — r h-uL^Eu <g> 
+(rii;jiLu — ru-jiLu — r^-uLu + ru;iiLji)Eu ® £yj 
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+ (rji-uLu - rn-aLji^Eji ® E u - E u ® E fl ) - (r ij;U L u - r K . K Iry) ® ~ ^» ® ^i)] 
+ ^ [(rii-jjLu — ru-jjLu + r^-uLji — Tji-uL^Eu ® E 1 ^- 

+ (ru-jjLu - ru.jjLii){Eu <g> Ejj - Ejj ® E u )] 

+ J] 4(r H; «L^ - ru-uL^Eu ® (2. 19) 



The first term of Eq. ( |2.19 ) is: 



(rji-stLu - r u . st Lji) E jt <g> E st 

i,j,s,t,l,s^t,j^l 

( r ji;stLu — Tu-stLji) Eji ® 

i,j,s,t,l;s^t,j,l;t^j,l;j^l 

i;fcZ-^ij — r ij;klLli) fyj ® E k [ + (V k i-ljL it — ru ; ljL ki ) E k i ® Eij 

i,j,k,l,jj^k,l;k^l 

+ {i~ji ; kiLu — ra-kiLji) Eji <g> + {ru-i k Lij — rij.^Lu) Ey ® Ei k ] 

+2 [{ru.jiLij - rij.jiLu) E^ ® + (r H;y Ly - r ij;lj L H ) Eij <g> . 

Therefore, we have 



{ r ji;stLu — rii. st Lji) Eji ® E st 



[ri2,Xi] - [r 2 i,L 2 

£ 

z,j,s,t,l;s^t,j,l;t^j,l;jjtl 

+ ^ [(ru ; kiLij — Tij-hiLu) Eij ® E M + (rji-ikLu — ru-^Lji) Eji ® 22, 

+ ( r ji-,kiLu — r u-kiLji) Eji ® E k i + (r u-i k Lij — r ij-i k Lu) Eij ® Ei k ] 
+ [(— J" ji;kkLu + r u- kk Lji + r ki -jiL ik — r ik -jiL ki )E kk ® Eji 

i,j,k,l;k^j,l,jytl 
+ ( r ji;kkLil — T u- kk Lji)Eji ® E kk \ 

+ [(^ii^ _ r ;z;Zj£zi + 2 {rij-iiLu - ru. tt Lij))Eu <g> J5y 



Ik 
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+ (th-jiLq - ni-jiLu + 2(-r ji . u L a + ru-uL^)) E u ® E jX 
+ (ru-uLij - Tij-uLu) Eij ®E u + 2 (/•/,:, //-,., - ny,jiL u ) E tj <g> ^ 
+ {rji-uLu - ru-uLji) Eji ®E u + 2 (/•/,:/., /-,., - n^La) E\j ® E l3 \ 
+2 (r li-jjLu — r u-jjLu — r muLy + r ij-uLji) Eu <g> Ej 3 

+4 ^(rii-uLu - r U; iiLii)Eii <g> E u 



i.i 



{Li, L 2 } by Eqs. (gTg ) and flCTQ . This finishes the proof of Eq. ( gig ). 
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